Spin injection through the depletion layer: a theory of spin-polarized p-n junctions 

and solar cells 
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A drift-difTusion model for spin-cliarge transport in spin-polarized p-n junctions is developed and 
solved numerically for a realistic set of material parameters based on GaAs. It is demonstrated that 
spin polarization can be injected through the depletion layer by both minority and majority carriers, 
making all-semiconductor devices such as spin-polarized solar cells and bipolar transistors feasible. 
Spin-polarized p-n junctions allow for spin-polarized current generation, spin amplification, voltage 
control of spin polarization, and a significant extension of spin diffusion range. 



Spintronics |jl| has played an important role in defin- 
ing novel applications that are either not feasible or inef- 
fective with traditional semiconductor electronics. Spin- 
tronic devices have found their niche in industries for 
magnetic read heads and nonvolatile memory cells. Here 
we propose and demonstrate a new scheme for spintron- 
ics, a spin-polarized p-n junction, which amplifies spin 
density, significantly extends the range of spin diffusion, 
electronically tailors spin polarization, and, in combi- 
nation with light as a spin-polarized solar sell, gener- 
ates spin-polarized currents with tunable spin polariza- 
tion. We prove these concepts by solving drift-diffusion 
equations for a realistic device model based on GaAs, 
which demonstrates that spin polarization can be in- 
jected through the depletion layer. Possibility of injecting 
spin polarization through a transistor is also discussed. 

The electrical injection of spin-polarized carriers 
within all-semiconductor structures (from a magnetic 
into a nonmagnetic semiconductor) was recently realized 
experimentally Q| (the scheme proposed in [^). Optical 
injection of spin-polarized carriers (both minority [^||^] 
and majority [^,^) has been known for some time. In 
addition, the relatively long spin diffusion lengths [^,0, 
coherent spin transport across semiconductor interfaces, 
a successful fabrication of a magnetic/nonmagnetic p-n 
junction Q based on the novel (Ga,Mn)As material 
and the recent demonstration of a gate- voltage control of 
magnetization in (In,Mn)As ||ic|] , make semiconductors 
promising materials for spintronic applications ||l^. 

In this paper we investigate the spin-charge trans- 
port in semiconductors under the conditions of inhomo- 
geneous bipolar doping (Flatte and Vignale |0 have re- 
cently made an interesting proposal for spin diodes and 
transistors in unipolar semiconductor heterostructures- 
a very different case from ours); we are not concerned 
with spin injection per se. Our model device is a spin- 
polarized p-n junction with spin polarization induced (ei- 
ther optically-in which case we get a spin-polarized solar 
cell-or electronically) to minority or majority carriers. 
By studying spin-charge transport numerically across the 
depletion layer, we observe novel phenomena, all result- 



ing from the fact that spin polarization is transferred 
(what we call injected) through the depletion layer. 

We first introduce a drift-diffusion model for spin- 
charge transport in an inhomogeneously doped semicon- 
ductor illuminated with circularly polarized light (and 
later solve the model for GaAs). In addition to the ap- 
proximations used in deriving the usual (unpolarized) 
equations [p^ , we assume that all the dopants are fully 
ionized, the carrier populations nondegenerate and vary- 
ing only in one (x) direction. Further, we assume the 
spin polarization to be carried solely by conduction elec- 
trons (that is, consider holes unpolarized), as is the case 
of III-V semiconductors like GaAs best candidates 
for photo-spintronics. 

The following parameters describe our model: accep- 
tor (donor) densities Na (Nd)', electron (hole) densities 
n (p), in equilibrium tiq (po); intrinsic carrier density 
rii (nf — riopo)] electron (hole) number current densi- 
ties Jn (Jp)', electron (hole) mobilities and diffusivities 
fin ifJ-p) and Dn (Dp); bipolar photoexcitation rate G; 
intrinsic generation-recombination rate constant w and 
spin-relaxation time Ti. Electron parameters n, J„, and 
G will carry spin index A (A — l,f for spin up and 
A = —1, i for spin down): -\- ni = n, J„| -|- J„x = J„, 
and G| + G| = G. We also define spin-related quantities: 
spin density s = rij — n^, spin polarization a — s/n, and 
spin photoexcitation rate G* = G| — Gj. 

The time evolution and spatial distribution of carrier 
and spin densities is described by three sets of equa- 
tions, (i) Poisson's equation cP(t)/dx'^ — —p/e, where 
is the electrostatic potential (related to the electric field 
E — —d(j)/dx in the x direction), e is the sample dielec- 
tric permittivity, and p = e{N£i — Na — n-\-p) is the local 
charge density with the elementary charge e. (ii) The 
linear response equations for the particle currents, J„a = 
—p,nn\E — Dn{dn\/dx) and Jp = -\-pppE — Dp{dp/dx). 
The mobilities and conductivities are connected via Ein- 
stein's relation eD = ksTp; ks is the Boltzmann con- 
stant, and T the absolute temperature, (iii) The conti- 
nuity equations 
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-w(nxp-nQpo/2) — KGa, (1) 



dt dx 
dp dJp 

— + -7- = -w(np - nopo) + G, 
dt dx 



2Ti 



(2) 



expressing the particle generation and recombination as 
well as spin relaxation. Equation sets (i-iii), together 
with appropriate boundary conditions drawn from the 
actual physical situation, fully determine the steady state 
{dnx/dt = dp/dt = 0) distribution of carrier and spin 
densities. 

In a homogeneous case and steady state Eqs. |l] and 
^ become wlrip — uoPq) = G, and wsp + s/Ti = G^. 
Note that the rate at which electron spin relaxes is not 
1/Ti, but rather wp + l/Ti, reflecting the disappear- 
ance of spin by recombination (spin is effectively trans- 
ferred to holes which then lose it). This is more pro- 
nounced in p-doped samples. Let us see if our equa- 
tions recover what is already known about spin po- 
larization in semiconductors. For polarization we get 
a = ao{l + uqPq / np) / {\ + 1/Tiwp), where ao = G" /G is 
the polarization at time of pair creation. Take a p-doped 
sample. Then p w p^, n <^ hq, and a — ao/{l + r/Ti), 
with electron lifetime t — l/wpo, expresses spin orien- 
tation [Q: if T ^ Ti, which is usually the case, a ~ ao 
and electron spins are effectively oriented throughout the 
sample. In an n-doped sample wp « G/uq, and the spin 
polarization is ao/(l + no/GTi). This is optical spin 
pumping jj]: spin polarization depends on the illumina- 
tion intensity G and is noticeable at G ^ uq/Ti. 

Our prototype model device is a 2 /xm long GaAs sam- 
ple, doped with Na = 3 x 10^'^ cm^'^ acceptors on the left 
and with Nd = 5 x 10^^ cm~'^ donors on the right along 
the X-axis (the doping profile is shown in Fig. The 
left surface of the sample, a; = 0, is illuminated by cir- 
cularly polarized light which creates electrons and holes 
with photodensities Sn = 6p = 3 x 10^^ cm~'^, and in- 
duces electronic spin polarization ao = 0.5 (the value 
given by the band-structure symmetry Q|) at the sur- 
face. Since holes lose their spin orientation faster than 
they gain it (also a band structure effect Q), their po- 
larization is effectively zero and need not be considered. 
The task is to find the steady-state distribution of elec- 
tron and hole densities n and p, as well as electron spin 
density s = n-^- — and polarization a = s/n, as a func- 
tion of X. 

We solve the equation sets (i-iii) numerically, in 
the steady-state regime, with the boundary conditions: 
p{0) = po(0) + Sp (where po{0) = 0.5{Na + {N^ + 
4n^)°-^)), n(0) = Po(0)/nf + Sn, and a(0) = 0.5, re- 
flecting bipolar photoexcitation and spin orientation pro- 
cesses at the illuminated surface (there is no illumina- 
tion in the bulk: G — 0); n{2) — rig (2) (where no (2) = 
O.biNo + iNl+Anjf p{2) = no{2)/nl anda(2) - 0, 
maintaining equilibrium at the right surface. For electro- 
static potential (p we use (/)(0) = — (fc^r/e) ln(po(0)/7T.i) 
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FIG. 1. Spin injection through the p-n junction. Electron 
and hole densities n and p are shown as a function of the dis- 
tance X from the illuminated surface (x = 0). The dashed 
line is the doping profile No — Na (scale not shown): it 
is p-type with Na = 3 x 10^^ cm~^ on the left and n-type 
with Nd = 5 X 10^^ cm~^ on the right; the transition region 
is between 0.25 and 0.55 fim. Also plotted is spin density 
s = n-f — ni and spin polarization a = n/s in the inset. The 
remarkable result that a extends well beyond the transition 
layer (and s is amplified) demonstrates both spin injection 
and spin density amplification. 



and(/)(2) = (fceT/e) ln(no(2))/ni)- V, where F is the ap- 
plied bias, conventionally defined with respect to the dark 
built-in value Q. The room-temperature parameters of 
GaAs used in our model are p^ : the intrinsic carrier 
density rii = 1.8 x 10^ cm~'^, electron and hole mobilities 
/i„ — lOfip — 4000 cm^ • V^^ • s^^, electron and hole dif- 
fusivities (from Einstein's relation) Z)„ — IQDp — 103.6 
cm^-s~^, and dielectric permittivity e — 13.1eo (where eo 
is the permittivity of free space). As for the generation- 
recombination rate w (assumed to be a constant inde- 
pendent of n and p), we take the value of (1/3) x 10~^ 
cm'^ • s~^. The electron lifetime in the p-region is then 
T « 1/(wNa) ~ 0.1 ns. Finally, spin-relaxation time Ti 
is taken to be 0.2 ns, a reasonable value for GaAs ||5|,|6|,|l5| . 

We first discuss the case oi V — (the applied reverse 
bias from an external battery cancels the forward volt- 
age due to photoexcitation), which is in Fig. 0. Elec- 
tronic density starts off with the value of n « (5n at 
a; = 0, decreases somewhat in the depletion layer, then 
rises by more than a decade to reach its equilibrium value 
of no « Nd in the n region. The hole density sharply 
decreases from approximately Na through the depletion 
layer until it becomes ~ nf/Nu, the equilibrium value 
on the n side. Spin density s is 0.5n at the illuminated 
surface, and essentially follows the spatial dependence of 
n through the depletion layer, but once in the n region 
it decays towards zero p6|. The surprising feature is the 



2 



increase of s through the depletion layer. In effect, the 
magnetization of the sample increases by more than an 
order of magnitude as a result of spin injection. The po- 
larization a (Fig. 1^, inset) starts at 0.5 at the illuminated 
surface, stays almost constant through the transition re- 
gion, then decreases to zero at the right boundary [ p^ . 
We checked that these results are robust against changes 
in Ti and r by up to two decades, as long as Ti ^ r (so 
that appreciable spin polarization can be induced in the 
p region in the first place). 

Spin polarization clearly survives the depletion layer. 
This is not an ordinary spin injection in which certain 
number of spin-polarized electrons tunnel through a con- 
tact and the spin density is equal on both sides. What we 
have, rather, is a spin pumping (leading to spin density 
amplification) mechanism. Indeed, in optical spin pump- 
ing circularly polarized light creates spin polarization 
of majority carriers (electrons in n region) by intensive 
illumination which generates spins at a faster rate than 
1 /Ti . Here we illuminate the p (not n) region, so we do 
not have optical spin pumping. The physics is the follow- 
ing: light induces spin polarization of minority carriers 
(electrons in the p region) through optical spin orienta- 
tion (see Ref. [|j), which diffuse towards the transition 
region where they are swiftly pushed by the built-in field 
into their native n region. What we have is spin pump- 
ing through the minority channel: spin-polarized minor- 
ity electrons bring spin into the n region faster than the 
spin there relaxes or diffuses away. 

In a sense, minority electrons play the role of circu- 
larly polarized light. The evidence for that is in Fig. |^. 
In the n region, spin imbalance is present even though 
the carrier populations are well relaxed in equilibrium. 
As a result, spin diffusion is controlled by the major- 
ity diffusivity constant _D„ |^ rather than the minority 
constant Dp, as would be the case of carrier diffusion. 
Spin decay beyond the depletion layer is therefore de- 
scribed by the single exponential s ^ exp(— x/Lg), where 
L" — (I?„ri)°-^ « 1.4 fim is the spin diffusion length of 
electrons in the n region. In the p-region spin diffusion 
length is = (DpTg)^'^; electrons diffuse bipolarly with 
hole diffusivity Dp for a time = tTi/(t + Ti) « 0.067 
ns, so that « 0.26 /im is of the order of the elec- 
tron diffusion length in the p-region. Since L"/Lf = 
{DyjDpf-^ X (1 -t- Ti/rf-^ (w 5 in our model), the ef- 
fective range of spin diffusion is extended far beyond the 
expected value. The depletion layer acts to extend the 
range of spin diffusion, with the effect most pronounced 
for Ti ;2> T, which is typically the case (especially at lower 
temperatures) . 

Our model scheme above constitutes a solar cell, as 
light illuminates the surface within the electron diffusion 
length from the depletion layer. As the electron diffu- 
sion length almost coincides with L^, electrons arrive at 
the transition region spin polarized. The built-in electric 
field then sweeps the electrons in the n-region (and holes 
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FIG. 2. Calculated spin polarization (q) profiles for differ- 
ent voltages. The upper figure is for the solar cell. The lower 
figure is for an injection of spin by a ferromagnetic electrode 
at j: = 2/im. In both cases a extends beyond the transition 
region, while the strong dependence of a on 1/ is seen only in 
the solar cell. The dashed line in both figures is the doping 
profile. 



back into p). Because of the spin amplification in the 
ri-region, the resulting electrical current is spin polarized 
and can be used for spintronic applications (in a combina- 
tion with ferromagnetic semiconductors or metals). Our 
spin-polarized solar cell has the usual I-V characteristics, 
with 1.03 V open-circuit voltage. 

How does spin polarization a(x) change with the ap- 
plied bias VI The calculated profiles for our model device 
are in Fig. |^ (top). There is a strong monotonic decrease 
of a with increasing V . To quantify this dependence 
we consider the change of the total spin accumulated in 
the cell (integral of the spin density s from a; = to 
2); the spin accumulates almost exclusively in the n re- 
gion (where electron density n is large). The result is in 
Fig. y. Total spin changes by almost 20 times when in- 
creasing V from -0.2 V to 1.2 V. By (loose) analogy with 
junction capacitance, we call this effect spin capacitance. 
The spin accumulation in the p region essentially follows 
the nonequilibrium density of electrons (charges) there, 
but in the n region, where spin diffusion length is much 
greater than carrier diffusion length, nonequilibrium spin 
accumulates to a much greater distance; spin capacitance 
is not trivially connected with nonequilibrium charges. 
Clearly, after switching off the light the spin is lost, so to 
"store" spins in the solar cell one needs to supply energy 
(the fact that Ti > r will lead to special transient effects 
related to spin recovery). The reason why a depends so 
strongly on V is that V changes the extent of the de- 
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FIG. 3. Total spin in the solar cell as a function of applied 
voltage. 

pletion layer: as V increases, the width of the depletion 
layer decreases [Q. Since we are illuminating the same 
point a; = in all cases, the amount of spin polarization 
that reaches the depletion layer decreases as the width of 
the layer decreases (so that the distance from the surface 
to the center of the layer increases). At large forward 
voltages (say, V = 1.2 V, in Fig. ||), the injected spin 
essentially decays within if, which is shorter than the 
distance from the illumination surface to the depiction 
layer. Another effect affecting the dependence of a on 
V is the value of the built-in field, which also decreases 
with increasing V. This electronic control of spin po- 
larization could be measured by observing luminescence 
of electrons in the n region (by, say, forcing them to re- 
combine with holes in heterostructures attached to the 
right surface). The above effects are not limited to op- 
tical spin injection. The only requirement is that there 
is an electronic spin and carrier imbalance (which can be 
also created electronically) in the p region. 

Finally, consider the case of spin polarization injected 
in the majority (here n) region (say, by ferromagnetic 
electrodes). There is no light illuminating the sample, 
and all the parameters and boundary conditions remain 
as above, except that at a; = all the charge densities 
are the equilibrium ones, and a{2) = 0.5 (a quite favor- 
able case). The calculated a{x) for different V are in 
Fig. H (bottom). The polarization, as before, is injected 
through the depletion layer (into minority electrons), but 
now V does not affect a{x) that much. Consider what 
would happen if an n region was now attached to the 
sample from the left, that is, we have an n-p-n transistor. 
The spin-polarization in the base (p) would be injected 
and amplified in the n-collector (by the spin-pumping 
from the minority channel we introduced above). We 
conclude that spin polarization can be injected all the 
way through a transistor, from emitter into collector. 
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